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Abstract 

The  cold  energy  of  liquefied  natural  gas  (LNG)  is  generally  wasted  when  the  LNG  is  extracted  for 
utilization.  This  paper  proposes  cryogenic  thermoelectric  generators  to  recover  this  cold  energy.  The  the¬ 
oretical  performance  of  the  generator  has  been  analyzed.  An  analytical  method  and  numerical  method  of 
calculation  of  the  optimum  parameters  of  the  generator  have  been  demonstrated. 
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1.  Introduction 

Thermoelectric  generators  are  driven  by  large  temperature  differences  to  convert  the  heat  flow 
through  them  into  electric  power  output.  Because  of  their  relatively  low  conversion  efficiency, 
their  applications  are  usually  limited  to  specific  situations  where  their  simple,  reliable  and 
compact  characteristics  present  advantages.  However,  one  exception  is  the  thermoelectric 
recovery  of  waste  heat  or  cold  when  the  large  temperature  difference  can  be  provided  without 
consideration  of  cost.  Consequently,  the  low  conversion  efficiency  is  not  a  serious  drawback. 
Natural  gas  is  known  as  a  clean  energy  source  commonly  used  as  a  domestic  and  industrial  fuel 
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for  combustion.  It  is  liquefied  to  transport  and  stored  at  atmospheric  pressure  and  at  a  tem¬ 
perature  of  about  -162  °C.  When  supplied  for  utilization,  the  liquefied  natural  gas  (LNG) 
should  be  vaporized.  Its  latent  heat  of  vaporization  and  any  sensible  heat  required  to  superheat 
the  vapour  are  termed  ‘cold  energy’,  which  is  generally  supplied  by  sea  water.  Cold  energy  is  a 
high  quality  energy  source  from  the  point  view  of  thermodynamics  [1].  If  thermoelectric  gen¬ 
erators  are  operated  during  the  vaporization  and  temperature  increasing  processes  of  the  LNG, 
with  the  two  ends  of  the  generator’s  arms  maintained  at  the  sea  water  temperature  and  the  LNG 
temperature,  respectively,  then  the  cold  energy  flow  though  the  arms  can  be  partly  converted 
into  electrical  power  output.  Previous  investigations  [2-5]  of  thermoelectric  generators  have  been 
focused  on  converting  waste  heat,  which  has  a  different  energy  source  and  operation  temperature 
range  from  that  for  converting  LNG  cold  energy.  For  thermoelectric  recovery  of  the  LNG  cold 
energy,  the  energy  source  is  the  heat  sink  with  temperature  well  below  ambient  temperature, 
which  provides  the  cold  energy,  namely,  absorbs  heat  from  the  ambient  environment.  The 
important  design  parameters  of  the  power  generators  are  the  efficiency,  power  output  and  the 
optimal  size  and  load  resistance.  In  the  present  paper,  we  firstly  give  the  analytic  solutions  with 
constant  transport  coefficients  and  compare  them  with  those  of  thermoelectric  generators  con¬ 
verting  waste  heat.  Then,  we  employ  a  numerical  method  with  real  temperature  dependent 
material  parameters.  The  calculation  results  of  the  two  methods  are  compared  and  analyzed  to 
give  a  profile  of  the  optimal  design  and  performance  of  thermoelectric  generators  converting 
LNG  cold  energy  to  electric  power  output. 


2.  Maximum  efficiency  and  power  output 


According  to  the  fundamental  theory  of  thermoelectricity  [6],  when  an  electric  current  density 
7  exists  in  addition  to  a  temperature  gradient  V?T  in  the  semiconductor  sample,  the  heat  flow 
density  (q)  and  the  heat  generated  per  unit  time  (q)  are  given  by: 

q=jTS-kVPT  (1) 

q  =  fp  ~  t(7  •  V?T)  +  V?(kV?T)  (2) 


where  S  is  the  Seebeck  coefficient,  k  the  thermal  conductivity,  p  the  electric  resistivity,  and  t  the 
Thomson  coefficient,  which  is  provided  by  the  first  Kelvin  relation  with  Seebeck  coefficient  as: 


dS 

X=Tdf 


(3) 


The  second  term  in  Eq.  (2)  is  the  Thomson  heat,  while  the  first  and  third  terms  are  the  Joule 
heat  and  heat  transport  by  thermal  conduction,  respectively. 

For  a  standard  thermoelectric  generator  as  shown  in  Fig.  1,  one  end  of  its  arms  remains  at  the 
temperature  Th  of  the  heat  source,  while  the  other  end  remains  at  the  temperature  Tc  of  the  heat 
sink.  The  current  is  assumed  to  flow  in  one  dimension,  as  along  the  arm  of  the  generator.  We 
assume  steady  state.  Therefore,  in  an  «-type  arm,  Eqs.  (1)  and  (2)  are  presented  as: 

q(x)  jT(x)S(x)  -  k(x)d^ 


(4) 
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Fig.  1.  Schematic  diagram  of  a  thermoelectric  generator. 


j2p(x)  -  T 


dS'(x) 


d2W 

dx  J 


(5) 


As  in  the  p-type  arm,  there  exist  similar  equations,  in  which  the  only  difference  is  that  the 
electrical  current  flows  in  the  opposite  direction  of  that  in  the  n  arm,  so 

I  =  jrAu  =  -y'pdp  (6) 

where  An,Ap  are  the  cross-sectional  areas  of  the  n  arm  and  p  arm,  respectively. 

For  calculation  simplicity,  the  material  parameters  ( S,p,k )  are  taken  as  constants,  or  if  they 
depend  upon  temperature,  the  values  are  taken  at  the  mean  temperature  of  the  materials.  With 
these  approximations,  the  solutions  to  Eqs.  (4)  and  (5)  are: 

T(x)=.Te  +  jAT  +  ^x(l-x)  (7) 


q(x)  =  jST(x )  -  k^j-  -  ^ pj2(l  -  2x)  (8) 

where  AT  =  —  Tc,  and  /  is  the  length  of  the  arms. 

The  temperature  distribution  along  the  arm  of  the  device  is  given  by  Eq.  (7).  We  can  then 
obtain  the  heat  flow  density  at  each  end  of  the  device  qx=0  and  qx-i  by  Eq.  (8).  From  the  energy 
balance,  we  find  that  the  power  output  of  the  generator  equals  the  difference  between  the  heat 
absorbed  from  the  hot  side  and  the  heat  supplied  to  the  cold  side.  In  using  the  LNG  cold  energy  to 
generate  electrical  power,  the  efficiency  of  the  generator  r\  is  defined  as  the  power  output  divided 
by  the  heat  supplied  to  the  cold  reservoir. 

_  On.x-iAn  +  qp,x=iAp  -  qa,x=0An  -  qp,x=0Ap  _  qB^i  +  qp.x  if)  -  qnx-0  -  qpx  0/i 
qnx-()An  +  qpx-()Ap  qn,x=o  +  Qp,x=oP 


where  fi  =  Ap/An. 
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Combined  with  Eqs.  (6)— (8),  Eq.  (9)  gives  the  efficiency  of  the  device  for  given  /  and  fi.  From 
=  0  and  p  =  0,  ij  can  be  maximized  as: 


_  AT(1  -  1/e) 
nmax  ~  Tc  +  Th/e 


(10) 


when 


AS  AT 
r(s  +  1) 


(11) 


(12) 


where  AS  =  Sp  -  Sn  (Sp  is  the  Seebeck  coefficient  of  the  p-type  material,  and  Sn  is  the  Seebeck 
coefficient  of  the  n-type  material  in  the  device),  r  is  the  internal  electrical  resistance  of  the  device, 
which  is  defined  as  r  =  ( pn/An  +  pp/Ap)l. 

e=y/l  +ZTm  (13) 

where  Tm  =  (Th  +  Tc)/2  is  the  mean  temperature,  the  figure  of  merit  Z  =  Afr,  and  K  is  the  thermal 
conductance  of  the  device,  defined  as 


K  =  kr 


An 

'T 


+  kp 


ap 

T 


To  have  ^max,  the  outside  load  resistance  should  be 


(14) 


The  maximizing  values  /?  in  Eq.  (12)  and  Rkad  in  Eq.  (14)  are  identical  to  those  for  waste  heat 
conversion  generators  [3],  However,  the  maximum  efficiency  for  the  later  is  r]maxheat  = 
Comparing  with  Eq.  (10)  indicates  that  the  efficiency  for  converting  LNG  cold  energy  should  be 
larger  than  that  for  converting  waste  heat,  which  is  reasonable  because  cold  energy  is  a  higher 
quality  energy  source  than  heat  from  the  point  of  view  of  thermodynamics.  Fig.  2  shows  the 
maximum  conversion  efficiency  of  these  two  types  of  generators.  We  can  also  deduce  from 
Eq.  (10)  that  if  the  dimensionless  figure  of  merit  ZTm  dose  not  change  much,  the  device  should 
have  better  performance  between  larger  temperature  differences  AT  and  in  a  lower  temperature 
region  since  the  efficiency  increases  with  the  decrease  of  Th  and  Tc.  As  e  — ►  oo  the  efficiency  ap¬ 
proaches  the  Carnot  efficiency. 

The  electrical  power  output  is  given  by: 

P  =  I(AS  ■  AT  -  Ir)  (15) 

At  the  maximum  efficiency,  /  can  be  eliminated  by  Eq.  (11),  and  by  replacing  r  by  its  definition, 
Eq.  (15)  yields: 
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Fig.  2.  Maximum  efficiency  of  thermoelectric  generators  vs.  temperature  difference  AT.  Tc  for  waste  heat  conversion 
and  Th  for  LNG  cold  energy  conversion  are  the  ambient  temperature  of  290  K.  The  figure  of  merit  Z  is  assumed  as  0.002 
K-1. 


/ASAT\2e  _An /ASAT\2  s 
V  e  +  1  /  r  ~  l  \  £  +  1  /  (Pn  +  Pp/P) 


(16) 


The  power  output  of  the  device  is  proportional  to  the  ratio  of  the  cross-sectional  area  to  the 
length  of  the  thermoelement,  which  suggests  that  the  shorter  the  thermoelectric  arm  the  larger  will 
be  the  power  output  per  unit  area,  consistent  with  what  is  shown  in  Fig.  3  in  Ref.  [4]  by  measuring 
several  thermoelectric  modules  of  different  sizes.  However,  for  real  thermoelectric  generators,  this 
is  not  always  the  case  because  of  the  unavoidable  thermal  resistance  between  the  arms  and  the 
heat  reservoirs  caused  by  the  conducting  strips  between  the  arms  and  the  electrical  insulating  layer 
at  each  end  in  addition  to  the  contact  resistance.  After  the  length  decreases  to  a  certain  value,  if 
the  arms  become  shorter,  then  the  temperature  difference  will  largely  occur  on  these  thermal 
resistances  rather  on  the  arms,  and  thus,  the  power  output  will  consequently  decrease  sharply. 


3.  Numerical  approach  with  temperature  dependent  material  properties 

When  thermoelectric  generators  convert  LNG  cold  energy  to  electric  power  output,  the 
working  temperatures  usually  range  from  ambient  temperature  to  below  ambient  by  more  than 
-150  °C,  and  the  material  parameters  vary  with  temperature  largely  in  this  temperature  region 
according  to  experimental  data  [7,8].  Therefore,  analysis  of  the  performance  of  the  device  by 
assuming  constant  material  parameters  must  bring  in  considerable  deviation  from  the  true  values. 
Whenever  the  transport  coefficients  vary  with  temperature,  numerical  procedures  must  be  used  to 
solve  Eqs.  (4),  (5)  and  (9)  [5],  By  algebraic  manipulations,  Eqs.  (4)  and  (5)  can  be  converted  into: 

=  [7 T{x)S(x)  -  q(x)\/k(x) 


(17) 
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dq 

dx 


=  p{x)f[\  +  Z{x)T{x)}- 


jS(x)q{x) 


A(.v) 


(18) 


where  x  =  x/l,  j  =  jl  and  q  =  ql.  For  a  given  j,  T(x)  and  q(x)  along  the  arm  can  be  obtained  from 
the  above  two  recursion  relations  with  the  boundary  conditions  T(x  =  0)  =  Tc  and  T(x  =  1 )  =  Th. 
The  efficiency  of  the  device  can  then  be  calculated  by  Eq.  (9)  as: 

_  ^I-I  +  ~  Vpx-qP  ^ 

On.x-o  +  qp-(-a[j 


Therefore,  the  maximum  efficiency  can  be  obtained  by  finding  the  optimal  values  of  j  and  [S. 
The  electric  power  output  can  be  calculated  by: 


PVax  =  (Pn  +  Vp)I 


where  r,  =  f*  dxp{(x),  i  =  n,  p. 


--mi: 


dx|5n(x)|- 


I-  J  dxl^p^ 


(20) 


4.  Discussions 


In  bulk  materials,  Bi  Sb  alloys  as  n-type  and  CsBi4Te6  as  p-type  material  have  been  reported  to 
have  the  highest  figure  of  merit  up  to  now  at  low  temperature  [9],  However,  a  complete  TE  device 
needs  both  a  p-type  and  an  n-type  version  of  a  material  to  operate.  Thermoelectric  modules 
composed  of  p-type  Bio.5Sb1.5Te3  and  n-type  Bi2Te3_ASeA  with  x  =  0.12-0.3  is  a  practical  material 
for  cold  energy  conversion  thermoelectric  generators.  In  the  analytic  method,  the  figure  of  merit  Z 
of  the  material  is  assumed  to  be  a  constant,  however,  for  real  materials,  it  may  vary  greatly  with 
decreasing  temperature,  dropping  from  0.003  K  1  at  300  K  to  0.0015  K  1  at  150  K  [7].  We  cal¬ 
culate  the  device  parameters  of  a  cryogenic  thermoelectric  generator  by  both  the  analytic  method 
with  the  value  of  the  material  transport  coefficients  ( S,p,k )  at  the  mean  temperature  and  the 
numerical  method  with  these  coefficients  being  temperature  dependent.  The  device  is  composed  of 
p-type  Bi2Te3-Sb2Te3  alloy  and  n-type  Bi2Te3-Bi2Se3  alloy  with  their  transport  coefficients 
polynomial  fitted  with  S(T),  \/p(T),  k(T)  curves  measured  by  Ref.  [8]  for  No.  1  materials. 


Sp  =  -48.27  +  1.342T  -  1.650  x  10-3T2 

Sn  =  -16.02  -0.9636T+  1.120  x  10-3T2 

l/pp=  1.314  x  104  —  99.15T  +  0.2925T2  —  3.082  x  10-4r3 

l/pn=  1.154  x  104  -  101.7T  +  0.3502T2  - 4.247  x  10-4T3 

kv  =  73.07  -  0.6057T  +  2.290  x  10~3:r2  -  3.068  x  10-6r3 

kn  =  45.05  -  0.247 IT  +  6.892  x  10-4r2  -  6.752  x  lQ-’T3 


fiV/K 
fiV/K 
(Qcm)  1 
(Qcm)-1 
mW  /  (cm  K) 
mW  /  (cm  K) 


The  calculations  results  are  listed  in  Table  1.  The  operating  temperature  range  is  from  Tc  =  130 
K  to  rh  =  290  K.  The  area  An  and  the  length  of  the  device  are  1  cm2  and  1  cm,  respectively. 

It  is  quite  interesting  to  note  that  although  the  material  parameters  change  greatly  with  tem¬ 
perature,  the  analytic  method  using  mean  temperature  values  gives  a  very  good  approximation 
relative  to  the  exact  numerical  solutions.  From  our  calculations,  we  find  that  the  material 
transport  coefficients  at  the  mean  temperature  are  close  to  the  average  values  for  the  operating 
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Table  1 

Values  of  relevant  parameters  of  the  thermoelectric  device  adjusted  to  maximize  the  LNG  cold  energy  conversion 
efficiency 


*7  max 

/(A) 

Qc  (W) 

-^load/r 

P 

Analytic  method 

0.096 

20.88 

6.43 

1.273 

0.84 

Numerical  solution 

0.090 

20.53 

6.64 

1.293 

0.86 

Tc  =  130  K,  Th  =  290  K,  No.l  materials  (p-type  Bi2Te3-Sb2Te3  alloy  and  n-type  Bi2Te3-Bi2Se3  alloy)  in  Ref.  [8]. 


temperature  range.  Although  thermoelectric  materials  have  a  variable  stoichiometry,  the  trans¬ 
port  coefficients  may  vary  wildly  with  different  fabrication  processes.  The  temperature  dependence 
of  the  transport  coefficients  of  the  No.l  materials  in  Ref.  [8]  is  typical,  and  therefore,  the  values  in 
Table  1  can  give  an  estimate  of  the  accuracy  of  the  analytic  method.  Ref.  [5]  also  shows  that  the 
analytic  method  represents  a  good  approximation,  although  the  material  properties  and  tem¬ 
perature  range  are  different  from  our  case.  It  is  good  news  for  engineers  that  once  the  temperature 
range  is  determined,  the  optimal  design  values  can  be  easily  calculated  with  the  materials  trans¬ 
port  coefficients  at  the  mean  temperature.  The  maximum  efficiency  given  by  the  analytic  method 
tends  to  be  higher  than  the  actual  value.  Nevertheless,  the  actual  conversion  efficiency  can  reach 
9%,  which  is  very  good  for  thermoelectric  generators  with  a  temperature  difference  of  only  160  K. 


5.  Conclusions 

The  performance  of  a  LNG  cold  energy  conversion  thermoelectric  generator  has  been  inves¬ 
tigated  analytically  and  numerically.  Results  show  that  the  efficiency  of  cold  energy  conversion  is 
better  than  that  of  waste  heat  conversion,  although  the  optimum  parameters  of  Ap/An  and 
external  resistance  i?ioad  have  the  same  calculation  equations  for  these  two  kinds  of  energy  sources. 
The  calculation  results  indicate  that  the  analytical  method  with  material  parameters  evaluated  at 
the  mean  temperature  provides  good  accuracy  for  optimum  parameters  engineering  design  for 
LNG  cold  energy  conversion  generators. 
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